Brownian motion and the hydrodynamic friction tensor 
for colloidal particles of arbitrary shape 
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We synthesize colloidal particles with various anisotropic shapes and track their orient at ionally 
resolved Brownian trajectories using confocal microscopy. An analysis of appropriate short-time 
correlation functions provides direct access to the hydrodynamic friction tensor of the particles 
revealing nontrivial couplings between the translational and rotational degrees of freedom. The 
results are consistent with calculations of the hydrodynamic friction tensor in the low-Reynolds- 
number regime for the experimentally determined particle shapes. 
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Spherical colloidal particles have served as a model 
system for investigating Brownian motion since the pi- 
oneering studies of A. Einstein P and M. J. Perrin p]. 
Such particles are characterized by a translational dif- 
fusion coefficient that is linked to the Stokes friction 
coefficient through the well-known Stokes-Einstein rela- 
tion [1 . Particles encountered in nature and industry 
usually have complex non-spherical shapes, and describ- 
ing their Brownian motion raises fundamental questions 
about how translational and rotational diffusion are cou- 
pled. However, in most studies, translational and rota- 
tional diffusion are considered separately, which is valid 
only for certain highly symmetrical particle shapes. 

Recently, model colloids with well-characterized but 
complex shapes have become available [3]-[6], which per- 
mits the quantitative study of the hydrodynamic cou- 
pling between translational and rotational diffusion for 
nontrivial particle shapes for the first time. 

In general, the dynamics of a colloidal particle sus- 
pended in a liquid is described by a Langevin equation 
that equates the Stokes friction forces and torques with 
random thermal forces and torques on a particle. For 
an arbitrary colloidal particle suspended in a liquid, the 
friction forces and torques are described by a symmetric 
second-rank hydrodynamic friction tensor 7^ [II |8] , which 
includes off-diagonal terms coupling the three transla- 
tional and three rotational degrees of freedom. In all, 1-L 
has 21 independent elements. For spherical particles, 1-L is 
diagonal, with two distinct entries corresponding to the 
inverse translational and rotational friction coefficients 
[9l[T0]. For rod-like particles, both the translational and 
rotational entries involve two different coefficients, corre- 
sponding to parallel and perpendicular particle orienta- 
tion, but 1-L remains diagonal, meaning that translation 
and rotation remain decoupled [HI [12] . 

For a general biaxial particle, H involves nonzero off- 
diagonal elements that couple translational and rota- 



tional motion. The corresponding Langevin equation 
involves intricate multiplicative noise terms due to this 
coupling, which makes a description of the Brownian dy- 
namics much more difficult. Although a first theoretical 
treatment dates back to P. F. Perrin [I3l[l5, it was not 
reconsidered until much later, and only by a few authors 
[71 [T5HT7] who never explicitly applied it to experiments 
for biaxial non-orthotropic particle shapes. 

In this Letter, we report experimental measurements 
and theoretical calculations of the hydrodynamic friction 
tensor for several anisotropic colloidal particles, including 
a general irregular biaxial shape with three fused spheres 
of different diameters. The particle shape and size are 
determined by confocal as well as scanning electron mi- 
croscopy (SEM). We track the Brownian trajectories of 
these anisotropic colloidal particles with full orientational 
resolution in real space by confocal microscopy. This 3D 
real-space technique allows for tracking the motion of 
arbitrarily complex colloidal particles, even in crowded 
environments. Based on the generalized Stokes-Einstein 
relation, we then propose a theoretical framework to ex- 
tract all independent hydrodynamic friction coefficients 
from the short-time limit of appropriate correlation func- 
tions. Our results are consistent with low- Reynolds- 
number hydrodynamic calculations of the friction ten- 
sor assuming stick boundary conditions of the solvent 
at the particle surface, where the experimentally deter- 
mined particle shapes are taken as an input. Since the 
full orientational resolution of the individual particle tra- 
jectories reveals the couplings between different degrees 
of freedom of Brownian motion, the information obtained 
by our analysis is much more basic and detailed than av- 
eraged quantities derived from dynamic light scattering 
[TF or sedimentation [19 experiments of biaxial colloidal 
particles. Our method can be used to analyze the Brow- 
nian dynamics of any rigid irregularly-shaped colloidal 
particles. 
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FIG. 1. (Color online) SEM images (gray) and images obtained by confocal microscopy (colored insets) of three different 
colloidal particles (regular trimer, regular tetramer, and irregular trimer) and the corresponding hydrodynamic friction tensors 
1-i determined in experiments and predicted for idealized particle shapes by hydrodynamic calculations using Hydrosub. The 
shown tensors are given in dimensions [7^^,^=1,2,3] = M-ni, ['H*^4 1'g] = [7^}=i',2',3] — M-^^^i ^ind [7^^,^=4,5,6] = ^im^. Since % 
depends on the center-of-mass position and orientation of the particles, the chosen coordinate systems are illustrated. 



For our experiments, we have prepared anisotropic col- 
loidal particles from fluorescently labeled, cross-linked 
PMMA (poly (methyl methacrylate)) spheres [20 using 
an emulsion-evaporation method [3^. The resulting clus- 
ter shapes are uniquely set by minimization of the second 
moment of the mass distribution [3 and can be confirmed 
by SEM. We have specifically chosen regular trimers and 
tetramers as well as an irregular trimer shown in Fig. [l] 
for their different symmetry properties. We idealized the 
particle cluster shapes as a composition of fused spheres 
and measured only the radii and relative distances be- 
tween the spheres. 

For the regular clusters, RITC (rhodamine B isothio- 
cyanate) labeled PMMA spheres 2.1 ±0.1 jitm (trimer) 
and 2.4 ±0.1 |LLm (tetramer) in diameter are employed. 
For the irregular trimer, 2.1 ±0.1 iitm and 1.3 ± 0.1 iitm 
RITC dyed spheres are combined with 1.7±0.1|LLm 
spheres labeled with NBD-MAEM (4-methylaminoethyl 
methacrylate- 7-nitrobenzo-2-oxa-l,3-diazol), which al- 
lows us to easily distinguish the different spheres with a 
confocal microscope. All sphere diameters are measured 
by static light scattering (SLS). The clusters are dis- 
persed in a TBAB (tetra butyl ammonium bromide) sat- 
urated '^78/22 (weight /weight) cyclohexyl bromide/ d5- 
decalin (CHB/decalin) mixture, which nearly matches 
the particles' density and index of refraction and has a 
dynamic (shear) viscosity 77 = 2.22 mPa s. The disper- 
sion is then put into rectangular glass capillaries (Vitro- 



tubes, 100 |LLm x 5 mm x 50 mm) and the ends are sealed 
with optical adhesive (Norland #81). 

The three-dimensional motion of the anisotropic par- 
ticles is observed using a Leica TCS SP5 confocal mi- 
croscope equipped with an argon laser (Ai = 488 nm 
and A2 = 543 nm) and an oil-immersion objective (Leica, 
63 X, 1.4 NA). The imaging speed is typically 70 stacks 
in z-direction per approximatively 0.8 s. All experiments 
are conducted at room temperate T = 294 K. The posi- 
tions of the individual spheres of each cluster are tracked 
using IDL routines (see Ref. [21 ). From these sphere po- 
sitions we calculate the center-of-mass positions, orienta- 
tions, and bond lengths of the clusters. For the regular 
trimer the bond length is 1.5 jitm, for the regular tetramer 
the centers of any two spheres are separated by 2.3 |Ltm, 
and for the irregular trimer the bond lengths are 2.2 jam 
between the big and the medium sphere as well as be- 
tween the big and the small sphere, and 1.7 \im between 
the medium and the small sphere. 

Apart from temperature T and solvent viscosity 77, the 
Brownian dynamics of a single rigid colloidal particle 
depends only on its shape and size, which enter in the 
6 X 6-dimensional symmetric hydrodynamic friction ten- 
sor 7/ [7[ E] . The latter relates the translational velocity 
V and the angular velocity Co of the particle to the hy- 
drodynamic drag force F and torque T that the particle 
experiences in the viscous solvent: K — —rjHS with 
K {F,f) and ? {v,uj). 




FIG. 2. (Color online) Time evolution of ten representative correlation functions Cij{t) for the irregular trimer determined in 
experiments (upper row) and predicted from a simulation (lower row) based on the idealized hydrodynamic friction tensor (see 
Fig.[l]). 



There are two possibilities for determining H for a 
given particle. It can either be obtained from its shape 
and size by a hydrodynamic calculation that involves 
solving the Stokes equation with stick boundary condi- 
tions for the solvent at the particle surface [8] , or it can be 
extracted from appropriate equilibrium short-time corre- 
lation functions. We have used the software Hydro sub 
[22l424| to follow the first route, where we used the exper- 
imentally determined particle shape, idealized by fused 
spheres, as input [25J. For a trimer and a tetramer of 
equal spheres as well as for an irregular trimer, results 
are shown in Fig. [l] [26]. For convenience, we have chosen 
the coordinate systems in such a way that the center of 
mass of a particle coincides with the origin of coordinates 
and the particle's planes of symmetry coincide with the 
coordinate planes, whenever this is possible. This choice 
of particle-fixed coordinate systems leads to a particu- 
larly simple structure of the hydrodynamic friction ten- 
sor with many vanishing non-diagonal elements j^. The 
remaining non- vanishing elements are highlighted in Fig. 

m 

The second route to access H is to measure the tra- 
jectory of the Brownian particle with full orient ational 
resolution, i. e., the combined knowledge of the center-of- 
mass position x{t) and the three mutually perpendicular 
normalized orientation vectors Ui{t) with i = 1,2,3 in 
dependence of time t. The key idea is now to consider a 



set of appropriate dynamical cross-correlation functions 

= (1) 

with z,j G {!,..., 6}, where ( • ) denotes a noise aver- 
age and the six-dimensional positional-orientational dis- 
placement vector X{t) = {Ax{t) ^ Au{t)) is defined by 
Ax{t) = x(t) - f (0) and Au(t) = \ Yl^.^i ^i(O) x ^i(^). 
where the latter is the appropriate expression for orien- 
tational displacements. The short-time limit of this set 
of cross-correlation functions gives access to the hydro- 
dynamic friction tensor 1-L via 

P=-lim^^, H=-^P"S (2) 

where V denotes the (generalized) diffusion tensor and 
/cb Boltzmann's constant. From this second route, based 
on the experimentally determined trajectories, we obtain 
the results presented in Fig. [T] The experimental results 
for 1-L are in good agreement with our hydrodynamic cal- 
culations; deviations are due to an idealization of the 
particle shape in the hydrodynamic calculations and due 
to the statistical error originating from the limited length 
of the measured trajectories. 

In reality, the particles are not compositions of per- 
fect spheres, but have rough surfaces and deformations 
near the overlap-areas of the spheres. While uncertainties 
in the size of the particles only lead to small deviations 
in the translational-translational elements 1-Lij=i^2,3 oc I 
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with the length scale /, these deviations are of greater 
relevance for the translational-rotational coupling ele- 
ments 'H^-^4'5'g (X and lead to large deviations in 
the rotational-rotational elements Hi,j=4,5,6 l^. Non- 
diagonal tensor elements, which should vanish by symme- 
try considerations (not highlighted in Fig.[T]), indeed have 
nonzero values in the experimentally determined friction 
tensor due to both irregularities in the actual particle 
shape and the statistical error because of the finite length 
of the trajectories. 

We finally address the full time dependence of the ba- 
sic cross-correlation functions Cij{t). Results for Cij{t) 
obtained from our measured trajectories of the irregular 
trimer particle are presented in Fig.[2j Figure |2|a) shows 
the translational mean-square displacements for the dif- 
ferent Cartesian components, which increase linearly for 
short times with the slope governed by the anisotropic 
short-time diffusion coefficients. At long times, there is 
another linear function in time, which is the same for all 
coordinates, since it is governed by the orient ationally- 
averaged long-time diffusion coefficient. Several cross- 
correlations between translational and rotational dis- 
placements are shown in Fig. [2|b). The absolute val- 
ues of these nontrivial correlations are initially increas- 
ing with time, but decorrelate for longer times. Finally, 
the rotational-rotational correlations shown in Fig. [2|c) 
are clearly positive and build up continuously as func- 
tions of time until they decay again and approach the 
same constant value for very long times. Here, we re- 
strict the presentation of correlation functions to this ir- 
regular particle, since it is the only particle with nonzero 
translational-rotational coupling elements and thus pro- 
vides the most nontrivial dynamics. 

An idealized time evolution of the correlation functions 
Cij{t) can be obtained from a given hydrodynamic fric- 
tion tensor Ti. With this tensor as an input, the Brow- 
nian motion of a colloidal particle can be simulated by 
solving its Langevin equation (see Ref. |T7]) numerically 
using a stochastic integrator of strong order 1.5 [27). Our 
results for the predicted correlation functions are pre- 
sented in Fig. [2] as well. A comparison with the experi- 
mental results for Cij{t) reveals again a good agreement 
with deviations resulting from choosing the idealized in- 
stead of the experimentally determined hydrodynamic 
friction tensor and from the statistical error. Note that 
the statistical error is obvious for the experimentally de- 
termined correlation functions, while it is extremely small 
for the simulated correlation functions, where trajectories 
with 10^ time steps have been calculated. 

In this work, we track the Brownian dynamics of indi- 
vidual colloidal particles with various anisotropic shapes 
and extract the hydrodynamic friction tensor from an 
analysis of appropriate short-time correlation functions. 
The framework of our analysis can in principle be applied 
to any rigid particle with an arbitrary shape and there- 
fore to a broad range of relevant suspensions. Using con- 



focal microscopy we obtain real-space 3D measurements 
of any complex particles, even in crowded environments. 

Future work should address the effect of aligning ex- 
ternal fields such as gravity or magnetic fields [28] that 
gain major importance in the context of directed self- 
assembly [29 . Moreover, the case of considerable colloid 
concentrations should be explored, where direct particle- 
particle interactions [30^ and solvent-fiow mediated hy- 
drodynamic interactions [31] will lead to even more intri- 
cate translational-rotational couplings. This would pro- 
vide a basis to understand the rheological properties of 
concentrated dispersions of irregularly shaped colloidal 
particles like clay [32] and asphalt [33] . 
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